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1 Introduction 



The Bern-Kosower method is described as a set of simple rules to obtain gluon scattering 
amplitudes at one loop , and it is known to improve the computational efficiency over the current 
Feynman diagram technique O . Those rules are derived from a string theory in the limit where the 
inverse string tension vanishes, as a consequence of the idea that a string world-sheet degenerates 
into a desired particle diagram at a singular point on the boundary of moduli space The 
integration over moduli space naturally covers all necessary Feynman diagrams appearing in field 
theory, and we hence have a compact master formula for particle scattering amplitudes. Thus the 
most conspicuous point in this formalism is that the diagram summation is already finished in the 
formula without introducing the loop integral and the Dirac trace for a given scattering 

!,!• This 

idea is also applied to graviton scattering {?]]. 

The discovery of the Bern-Kosower rules has also stimulated investigations for a new mathe- 
matical structure of quantum field theory; how to reflect the string-like structure into field theory 
as such. The first rederivation of the Bern-Kosower rules was accomplished by Strassler in Ref. |8| , 
where the background field method, the proper time method and the path integral method for a 
first quantized + 1 dimensional field theory (world- line formulation) are well combined || . There 



are many other fruitful examples along this stream |LC]] (see Section 1 of Ref. [11] for updated 
references), and these examples are the strong incentives to study the world- line formalism from 



the theoretical point of view, especially from the viewpoint of its higher loop extensions [12] . The 
present paper also discusses two-loop Yang-Mills theory in the world-line formalism with the aim 
to develop techniques for higher loops. 



This paper is a continuation of the previous work [11], where the effective action of Yang-Mills 
theory at the two-loop order is derived based on the world-line formalism; also developed there is a 
certain technique, which generates multiloop generalizations of the one-loop trace-log (determinant) 
formula. However, these arguments are still inside the shell of formal arguments, since we are left 
with the problem of how to deal with the multi-integrals of world-line moduli parameters. The 
moduli integrals of a higher genus world-sheet are too complicated to perform, while one might 
naturally expect that this situation would be improved in the field theory limit. Although we, of 
course, have an option to computerize these complicated integrals, there are still difficulties, for 
example in three loop QED integrals fl3| ]. 

It is certainly valuable to analyze two-loop integrals of Yang-Mills theory in the world-line 
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formalism in particular if many outer particles or an Euler-Heisenberg type constant field are 
involved. It might also hint to the world-line moduli integrations at higher loop orders. 

In this paper, after a detailed derivation of the two-loop Euler- Heisenberg action with gluon 
loops in a pseudo-abelian gauge field background, we shall present some technical issues of how to 
deal with the world-line moduli integrals in the gluon effective action at the second order of the 
Taylor expansion in terms of external background fields. This analysis is also essential to examine 
the divergence structure related to a wave function and gauge fixing parameter renormalization. 
We only discuss the gluon loop part, since the ghost loop part is rather simple and can be dealt 
with in the same way as the gluon loop case. 

When we perform the integrals, we insert a mass parameter in order to regularize divergences. 
Generally speaking, massive propagators in the Feynman rule method are difficult to integrate in an 
analytic way. Contrastingly in our formalism, we shall go through the entire procedure analytically, 
and all the results will be written in hypergeometric functions. This is certainly an intriguing point 
of this paper. 

This paper is organized as follows. In Section ||, we write down our starting formulae for the 
gluon loop effective action at two loops. We slightly modify a few notations from the previous pre- 
sentation [ 11 1 through the path and proper time inversions presented in Appendix |A|. In Section |3|, 



a derivation of the one-loop /^-function coefficient serves as an example for how calculations in our 
formalism can be simplified by specializing to the pseudo-abelian case. In Section we apply the 
pseudo-abelian technique to the calculation of the gluon effective action presented in Section |2[ 
Here, we only perform the world-line path integral parts. This Section is a completion of the parts 
outlined in the previous paper |f["l] ], and the details of the computation are contained in Appendix [B]. 
In Section [|, we further study how to integrate the world-line moduli integral parts, which are the 
final integrations to obtain a fully integrated form of Euler-Heisenberg type action. For simplicity, 
we only consider the gluon kinetic term, through the Taylor expansions concerning the external 
field strength. The Taylor coefficients are the functions of world-line moduli parameters, and we 
show that these coefficients can be integrated; the details are in Appendices ^ and [D]. Appendix |E| 
is the Feynman diagram analysis to be compared with our results. 

2 Two-loop effective action 

Let us first review in brief the world-line representation of the two-loop effective action in 
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Yang-Mills theory [11]. In this paper we only discuss the gluon loop part. It is given by 



V[A]=h[A\+h[A} 



(2) 

where (Ii = T\ + I^, 1% = I3 in the previous paper) 

00 S 00 w{T 3 )=x{r a ) 



and 



h[A) = ~ fdS fdr a fdT 3 f[Vx] s I [Vw\ Tz 







w(0)=x(0) 



K(0) - ±„(0)) w p (T 3 ) W^ p]u [x; S, r a ,0] 
+ x u (0)w p (T 3 ) W^[x;S,r a ,0}] W%[w;T 3 ,0] 



00 

h[A] = IJdTx dT 2 f[V Xl ] Tl j>[Vx 2 ] T2 <5 4 (x^(0) - 4(0)) 


x Tr c \\ a W p ,[ Xl ; T x , 0]1 Tr c [A a ^[x 2 ; T 2 , 0] 



with the following compact notations: 



[£>x]t F[x] = / Dxe F[x] for any functional -F[x] 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



( T 3 

W%[w; T 3 , 0] = P exp If drM T [w] 



(2.5) 



and 



W^ pu [x;S,r a ,0]=Tr c 



\ a P exp J y drM T [x] 1 A e P exp J f drM T [x] 



pv 



(M T [x))t = 2i 



F^Mr'))-^ v A c p (x(r'))-d T/ x^r') 



(A 



c\ab 



(2.6) 



(2.7) 



Here we have slightly changed the notations used in the previous paper [pTlJ] ; (i) the previous 
definition of M is associated with D = d — iA, while the present M T is associated with D = d + iA, 
(ii) the path ordering directions are modified to be the standard one, and some related formulae are 
listed in Appendix [A]. As in the previous paper, we always use Euclidean space-time conventions. 
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(a) 



(b) 



Figure 1: (a) The loop type parametrization. (b) The symmetric parametrization. 



It is convenient to have another representation for I\ [A] , based on the symmetric parametriza- 
tion in Figure 1(b), which treats the individual gluon lines in a more equal way, thus allowing a 
greater class of transformations by inverting and relabelling the gluon paths and leading to signif- 
icant simplifications in the concrete calculations. The expression (|2.2| ) is based on the loop type 



parametrization in Figure 1(a). The transformation rules between (a) and (b) are known |12| , 14 1, 
and we thus have the following symmetric representation for I\ [A] : 





r 3 x k (T k )=y 2 , 




n j[v*kH 




lk=1 x k (0)= yi 



(x 3/ ,(0) -xi M (0))x 3p (r 3 ) W% ap]v [xi Vi;T 2 , 0,11,0] 
+ x lu (0)x 3p (T 3 ) W'^[xJ 1 ) x 1 ;T 2} j ri > 0] ] W%[x 3 ;T 3 ,0] , 



where 



and 



S T a , 



w ae 



( T 2 \ ( Ti 

[x2~Vi;T 2 ,0,Ti,0] = Tr c A a P exp J JdrM T [x^ 1 ] > A e P exp J JdrM T [x x ] 



(2.8) 



(2- 



(2.10) 



Note that the transition from loop type to symmetric parametrization requires both splitting the 
loop path into two parts and inverting one of them (which we denote by x^ 1 ) to achieve all three 
paths to start at y\ and to end at y 2 - In contrast to the naive expectation, this suggests that 
in a general background, one may not just write down the product of three propagators starting 
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and ending at identical points to represent a loop with inserted propagator. See Appendix [A| for a 
detailed definition of the notation x^ 1 and some comments on path inversion 

After finishing the next section, we shall discuss the world-line path integrals for in Section |||, 
and the integrals of the world-line moduli parts (proper times) T)% in Section |5[ 



3 The pseudo-abelian case 

For the rest of this paper, we confine ourselves to the pseudo-abelian su(2) with constant field 
strength. Thus we assume 

A a Jx) = Ap(x)n a with n a n a = 1 and constant n , (3.1) 

i.e. the color dependence of the non-abelian gauge fields is supposed to be factored out in form of a 
constant unit vector in color space. Within these settings calculations are simplified considerably, 
though non-abelian results can still be reproduced as shall be seen below. As an example and 
to introduce some notations, we here show a brief sketch of how our formalism works within the 
calculation of one-loop /^-function coefficients. 



The assumption (3.1) leads to similar decompositions for the field strength 

F%{x)=F» v {x)n a with Tp V {x) = dpA u {x) - d v Ap{x) (3.2) 
and the matrix M T [x] 

M T [x] = M T [x] (n c A c ) = Mr [x] <g> T_ , (3.3) 
where A^T-fx] is the Lorentz matrix defined by 



M T [x] = 2i F{x) - - A p (x) x p 1 L (3.4) 



and 1l denotes the unit Lorentz matrix (in the Euclidean space). 

So far we have not used our additional assumption of a constant field strength, nor have we 
fixed the gauge for the external gauge fields A^\. If we take into account the constancy of the field 
strength, we may choose 



Ap,{x) X^J- ,yp , (3-5) 



henceforth expecting A4 T [x] to be of the form 



M T [x] = 2i F - - x a T up x p 1 L , (3.6) 



rather than (3.4). In addition, it is convenient to define the integrated matrices (omitting the index 
r), and we can simply write 

T T 

[ M T [x]dT = M[x}®T- with M[x] = J M T [x]dT . (3.7) 



It is the benefit of confining ourselves to the pseudo-abelian case, that the M T [;c] matrices for 
different values of the parameter r become commuting quantities 

[M T [x], M T t[x]\ =0 . (3.8) 

Thus we are allowed to drop path ordering from all of our expressions. This leads to the following 
decomposition 

P exp {M. <S> 71} = 1l®T + sinh{.M}<g>7L + cosh {M} ® T + , (3.9) 
in terms of the su(2) matrices 

71 = n c A c , T+ = (T_) 2 , 2=1 C -T + , where l c = diag(l, 1, 1) . (3.10) 
Now, using the properties 

Tr c 71 = 0, C A = Tr c T + = 2, Tr c J=l (3.11) 
the one- loop effective action for the gluon loop is calculated as follows [jl5|] : 

oo T 





oo 



2 J T € 
o 

oo 

1 f dT 



r^[A] = - 1 -J d fe-^ T f[Vx) T (Pe^jM T [x)dr)2 

o 

' j>[Vx\ T (D + C A Tr L {c.os\xM)) 

4 / ^fe~ m2T f[Vx] T (D + C A Tr L (e M )) + {? -» , (3.12) 
o 

where we have introduced the gluon mass term e~ m2T for regularization. The second contribution 
with T replaced by — J- counts for a factor of two, thus with the one-loop path integral normaliza- 
tion 

T 

j>Vx expj-i Jdr [x 2 + 2ixTx\ \ = (4^T)- D / 2 det^ 1/2 ( ) Jd D x (3.13) 

o 
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we are led to 



t1— loop 



[A] 



2(4tt)°/ 2 



drT _!_ D/2e _ m 2 T 



D + C A TV L (e 2 ^ J )det~ 



/ sm y - j \ 



d D x . (3.14) 



For now we are interested in the two-point function only, i.e. in the second functional derivative 
of the effective action with respect to A^. To this end, we only need the second order term of an 



expansion of ( 3.14 ) in terms of T . Omitting constant and higher order terms we find 



-il— loop 
G 



[A] 



oo 

|^(g-2) jdTT^f'e-^ Ji%^ + ... . 



2(4ir) D / 2 

U 

Setting D = 4 — 2e and performing the T integration leads to the following pole structure in s: 



(3.15) 



r 



1— loop 
G 



A] 



(it) I d ° X0 + o(£0) + • • • ' (3 - 16) 



where we have revived the gauge coupling g and where go is the dimensionless coupling constant 
defined by g = gofJ- e - Finally we calculate the functional derivative and transform into momentum 
space: using 



d D XQ 



1 



we read off 



n<# = (|) ( ^ - k,K ) + 0(e°) • 



(3.17) 



(3.18) 



Similarly, the one-loop contribution from a ghost loop can be calculated: As can be deduced 
from Ref. [11|, the ghost one- loop action is given by changing the overall normalization in ( 3.12| ) 
from 1/2 to —1, and only employing the Lorentz scalar term in ( |2.7| ); i.e., define —iA^x^X c = M T 
instead of using M T . The corresponding pseudo-abelian quantity M. is defined analogically to the 
gluon loop case (q.v. Eqs. fl3.6| ) and (3.7)). Thus we find the one-loop ghost action 

T 



■pi— loop 
L FP 



P [A] 



dT 
— < 
T 



dT 



-m 2 T 



-m 2 T 



1 f dT 



<£[Vx] T (Pexp J M T [x] drj 
o 

j>[Vx] T (l + CUcoshA^) 

T j[Vx] T {l + C A e M ) + {? 



(3.19) 



Again taking into account the —T term by a factor of two and using (3.13), we arrive at 



-il— loop 
FP 



[A] 



(4tt) d / 2 



dTT -l-D/2 e -m*T 



1 + CUdet L 



i/ 2 ( sin FT \ 
TT ) 



d D x • 



(3.20) 



Expanding this expression in the same way as done in (3.15), we derive 

oo 



and hence 



ab 



9 2 C A 8 a 



3e 



(3.21) 
(3.22) 

(3.23) 



Gathering Eqs. ( 3.1S| ) and ( |3.23| ), the correct (one-loop) /3-function coefficient 11/3 is reproduced. 



4 Two-loop Euler-Heisenberg formulas 

Now, let us consider the extension of the above calculations to the two-loop case. We deal with 



the symmetric representations ( 2^ ) and (|2.q ) . In the two- loop case , as understood from ([2/i]) , we 
have to keep in mind that the sign of the xTx term changes due to the r derivative, if we invert 
the path. For example in ( |2.10| ), one should notice (see also Appendix [A]) that 



(M T [x^\ 



ab 



2i 
2i 



^0*2 V)) ~ ? V^(*2 V)) • d T ^ p \r') 



t'=T2—t 



(A 
(A 



c\ab 



c\ab 



(M T2 _ T [x 2 ] 



6a 



(4.1) 



Reflecting this fact, it is rather convenient to introduce the signature index k (= ±1,0) on the 
Lorentz matrix M : 



M 



(k) def. 



dr2i 



. K X\ ta J : 7 p Xkp \lj 



(4.2) 



where the k stands for the line labels 1,2 and 3. With this matrix notation (f4.2j), the general form 
for the color matrix part of the action (^^) is written in the form 



W"§ a 8[xz Vi; T 2 , 0, T X , 0] W£[x 3 ;T 3 , 0] 



Trc 



A a exp [m { 2 ] ® 1- } g A e exp {m[ +) (g> T_ } exp {m ( 3 +) ®T_} 
After using the expansion (cf. Eq. ( |3.9[ )) 

exp [M ( k K) ®tJ} = 1 l ®1 + sinh [m[ k) }®T_ + cosh [m[ k) } ® T+ 



CO 



(4.3) 



(4.4) 
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we perform the color traces applying the following formulae: 



Tr c (\ a T±\ e T±)l ea = 2 
Tr c (\ a T± X e l) Tl a = 2 
Tr c (A a J \ e T±)Tl a = ±2, 

where the 3rd formula follows from the 2nd one with the properties 

(T ± ) T = ±T ± , 1 T =1. 

For any other combinations of A, B and C chosen out of {I, 71, 7+}, the following 
applies 

Tr c (\ a A \ e B) C ea = . 



cru 



(TV 



Thus the quantity Q4.3p is calculated as follows: 

T 2 , 0, T u 0] tOx 3 ; T 3 , 0] 
2 5 a/3 cosh {M - *} 5 cosh {.M 3 +) } ^+ sinhj.M^} a sinh{^ +) } 
+ 2 5 7 5 cosh cosh {M { 3 +) }^- smh{M[ +) } ^sinh^M^} 

+ 2<5 (7 Jcoshf7Wi +) } coshfTW^H + smh\M{ +) } sinh{^l _) } 

= ex P {A4f}^ex P {>(^} 75 ex P {M +) }^ 
+ exp{^W} a/3 exp{^f} 75 exp{-M +) }^ 
+ exp{^i +) } expj^-H ex P {.Mf } + ^ -T) . 

Now, as sketched in Ref. [11], performing the trivial integral (the 1st term) in Eq.(|4,2|) 

T k 

e^{fdr{-\xl)}^{±M[ K) }^ 
o 

1 Tfe 

= exp {±2i \k\ T k F}^ exp|-- Jdr[x 2 k + 2i (±k) x ka T ap x kp ]j , 

o 

and introducing the quantities 



3 T * 



5(«i.« 3 ,«8) = _ £ \dT\x\ + 2iK k x ka T (jp x k( \ , (« o = ±l,0) 



fc=i 
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we have the following formula: 



3 fc 

c p{"iE / dTxl}w^[x^\x 1 ;T 2 ,0,T l ,0}W^[x 3 -,T 3 ,0} 
k=l o 

= 5 a p exp{2iT 2 J 7 } exp{2iT 3 J 7 }^ 

y J a/3 I. J ua 

Sou expteiTif] exp(2iT 2 ^} e s<+.-.°) + _, 

L J a/3 I J -yd 



s s(°.-.+) 

5 (+,o,-) 



+ <$«„, exp^t'T!^ _exp42tT 2 .FK e °^""" + -»• - :F) . (4.12) 

It is worth noticing here that the interaction terms x^Tx^ defined on the three different lines possess 
different k values; this fact is related to the su(2) structure e abc . In the following, we calculate I\ 
and I 2 separately, since these two quantities involve different world-line topology. 

4.1 The h[A} part 

Applying the formula ( 4.12j ) to Eq. ( |2 . 8| ) , we obtain the following expressions (For the con- 
venience of presentation, we split h[A] into two quantities depending on whether x\x 3 or x 3 x 3 
correlations.): 

h[A] = T 1 [A]+T 2 [A] , 



where 



oo 

Ti[A] = ~JdT 1( lT 2 dn [d D yi Jd D y 2 



n 

.k=l 



e 2iT 2 F ( e KT 3 F_ lLTlL ( e 2iT 3 ^ 



Xk(T k )=y2 ■ 

T>Xk 

«(o,-,+) 



4(0) xg(T 3 ) 



hp 



+ 
+ 



( 



2iT 2 F ( '-HTxT 



l L Tr L (e 



-2,(T 1+ T 3 )^_ lLTri(e - 



2i(T!+T 3 )^ 



jg(+,-,o) 



c(+,o,-) 



(4.13) 



+ {F - -^) , (4.14) 



oo 




" 3 




Jdndndn Jd D yi j 




n 









.fe=i 


ccfe(0)=2/i 



^(O) 4(T 3 



-2isin(2^(T 2 + T 3 )) - e»Ci-^ + ^mr^^r 
+ [-2i sin(2^(Ti + T 2 )) - e^-W + gM^^M^ 
+ [-2cos(2^(T 1 + T 3 )) + e 2 ^-^ + i L Tr L e 2l ( Tl+T ^ 



<;(o,-,+) 



MP 



«(+,-,o) 



,S(+,o,-) 



(4.15) 



11 



Then we perform the path integrals of the form 





- 3 


Xk( T k)=y2 -, 


fd D yi Jd D y 2 


n 






_fc=i 





(4.16) 



and this yields 



where 



< i£(r) x v b (r') > (Kl , K2A 3)=-AA (Kl ' K2 ' K3) ^^^(r,r'; Kl , K2 , K3 ) , 



(4.17) 



M (^s) = {A7T) -D det -V2 K ^ cQt ^ ^ Tfc ^/ 2 det -i/2 , (4.18) 



and 



0#(r, t 7 ; «i, K 2 , «s) = -<5a6G^(r, r') 



fc=i 



2iK, a Tr _y \s 



1 1 



(4.19) 



with 



G^(r,r') 



VG , |(r,r') = <5 ( 
J-f ^ 



IT — T 



, 



T — T 



e-^^^ G -^') +m a ^ T G^(T,T') - i-) 



(«« = 0) 
(«« / 0) . 

(4.20) 



Inserting each value of (4-17) at (t, t') = (0, T3) into Eqs. (4.14) and (4.15), we therefore obtain 
(the details are presented in Appendix ||) 



00 , 

h[A) = ~ (4tt)- d Jdndndn det 1 / 2 ( 



TV L | 



2sin^"Ti cos2^(ri + 2T 2 ) -2smf(T 1 + T 2 ) coe.F(2Ti + 3T 2 ) 



Ajr sin TT 2 

+ { 1 - 2 cos aF(Ti + T 2 ) } sin ^"T 2 cos .F(Ti - T 2 ) 
J* 7 



+ 



4 sin .FTi sin .FT 2 sin 2TiT x + T 2 ) - 2 sin JFTi cos F(2Ti + 3T 2 



2 sin TT 2 cos ^(Ti - 2T 2 ) - sin F(T X + T 2 ) cos ZF(T a - T 2 



+ TV L | 



Ajr sin TT 2 



sin J"(Ti + T 2 ) cos^"(2Ti + T 2 ) - sin^Ti cos ZF(Ti + T 2 ) 



1 In Eq. (4.19), one may replace G^(r, r') — > G^ v (t, t') — G^(r, 0) - G£„(0, r') as seen in however our final 



results do not change because we only need the derivatives drdr'QI 
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+ 



T r 



r 2 n 



3 sin TT X cos T(2T X + T 2 ) + cos 2TT X sin T{T X + T 2 ) • Tr L ^cos 2FT 2 ^j 



Ajr sin JFT 2 



sin TT 2 cos J-"(Ti — T2) + cos TT 2 sin JF(Ti + T 2 



— sin TT\ cos 2.7 r T 2 + — — sin TT\ cos • TrH cos 2J T (Ti + T 2 ) 

r J \ j 

+ 5(T 2 )2(1-D)Tt l (cos2TT 1 \ + 5(T 3 )Tr L (cos2F(T 1 - T 2 )\ 



-S(T 3 )Tt L lcoB2jrr 1 j •Tr L ^cos2^T 2 j \ \d D x , 



where 



Ayr = sin TT\ sin-FT^ + TT$ sinjF(Ti + T 2 ) . 



(4.21) 
(4.22) 



4.2 The J 2 [A] part 



The computation of the other quantity I 2 [A] is similar to the above calculations, however the 
topology of the world-line diagram is different in this case. Let us start with the following expression. 
First, Eq. fl2.3|) with ( |3.7D inserted becomes 



00 

h[A] = ~ / dT^T 2 i[V Xl ] Tl j>[Vx 2 ] T2 5(xi(0) - x 2 (0)) 

A a exp [M { 2 +) ® 71 } 





xTr c 



A a expfA^S +) ®T_} 
With the expansion (4.4) and the properties 

Tr c (X a T + ) = Tv c (A a X) = 



UfJ, 



(4.23) 



(4.24) 



Tr c (A a 71) = n b Tr c (A a A b ) = n fe 2 5 ab = 2 n a , 



we have the formula 



Tr 



c 



(4.25) 



A a exp{^ K) (8)r_} = 2n a sinh{7WP}^ . (4.26) 

Remembering the relation n a n a = 1, we then derive from ( |4.23| ) 

00 

M = JWl jlV Xl]Tl f\V X2]T , S M 0) - „(0)) tr t [sinh {M^} stah^Wf 


00 

1 



(4.27) 
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where 



M K1 ' K2) = U D yi / d D y 2 



- 2 

n 

.fc=i 



Xk(T k )=y2 



Vx k 



5 D (yi-y 2 )e s ^^ 



with 



2 ±k 
*=io 

The normalizations J\f( K1 ' K2 ") satisfy the following properties: 

^+,+) =AA (+,-) =AA (-,+) =AA (-,-) 5 



(4.28) 



(4.29) 



(4.30) 



since, in the present case, the inversions of paths (i.e., the changes of re's signs ) do not change the 
value of the path integral ( |4.2Sj ): note that all the initial and ending points of two closed loops are 
identical. The (J- — > — J-) terms in ([4.27|) lead again to a factor of two, and thus we have 



(4.31) 



oo 

h [A] = I JdTidTzN^ Tr L (e 2t ^ +T ^ - 



The quantity (4.28) can be evaluated as follows. Recalling the relation 



8 (2/1 — 1/2) = Jim 



. D 

1 \T 



exp 



1 



T 3 ^o V47rT 3 

we convert the 5-function to the following path integral form: 



213(33) =3/2 



(5 {yi - y 2 ) = lim 

T 3 ^0 



1 



r 3 



£>x 3 exp[-- / x 3 (t)cIt 



and this leads to the relations 



= Bm^) = (4vr)~ fl detZ 1/2 ^ - ^ — ^ ^ 



T 3 ^0 

We therefore have the expression 

1 



^2 [X] = - (4vr) - D /" tffi (iT 2 det~ 



^ 2 



1/2 / sin TT\ sin .7^2 
^2 



V 



d xq 



Tr L cos 2J^(Ti + T 2 ) - Tr L cos 2T{T X - T 2 ) ld D x 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



It is interesting that the normalization J\f( Kl < K2 ) can still be obtained as a singular limit (the T3 —* 
limit) from Eq. ( 4.18; ), although /2 is not a singular part of ii (cf. ( |4.21| )) JlT| . 
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5 The world-line moduli integrals 

In this section, we study the divergence structure of the Euler-Heisenberg-type action derived in 
the previous sections. In the world-line formalism, divergences come out explicitly after performing 
the proper time integrals (the T a integrals; a = 1,2,3 in the present case), and hence we have 
to examine how to perform these integrals before discussing renormalizations. However, generally 
speaking, multi-integrations are difficult to perform in an analytic way, and hence we here consider 
the simplest case corresponding to gluon two-point function parts. 

Let us consider the Taylor expansions of T[A] concerning T (omitting constant terms) in the 
same way as done in the one-loop case. Expanding I\ and I 2 up to the second order of we 



extract the following quantities from Eqs. ( 4.21 ) and ( 4.35 ): 

„4 4e p no 

h [A] = (4 y 4 _ 2e JdPxQ T^ Vi \ 30Ci + 42C 2 + 48C 3 + —C 4 - 2C 5 - WCq 

57 17 
- (17d + 29C 2 + —C 3 + 32C 4 - — Cq) e 



- (2d + 2C 2 + 5C 3 + 6C4 - ^C 6 ) e 2 } + • • • , (5.1) 
h[A\ = ■^^ 2l Jd D x T^(-A)C 5 + ---, (5.2) 



with 



d = J dT\dT 2 dT 3 A e ~ 4 T^T 2 e ~ m ^ Tl+T2+T ^ 


00 _ 
C 2 = J dT\dT 2 dT 3 A £_4 TfT 2 e~ ( Tl+T2+T3 ^ 


C 3 = J dT 1 dT 2 dT 3 A £ ~ i T?T 2 T 3 e- m2 ( Tl+T2+T ^ 


00 „ 
C 4 = / dT\dT 2 dT 3 A e ~ 4 T^T 2 2 T 3 e ~ m ( T i+ T 2+ T 3) 


C 5 = J dT\dT 2 dT 3 A e_4 TfT$ S(T 3 ) e ~ m ( Tl+T2+T ^ 



(5.3) 





00 



C 6 = JdT 1 dT 2 dT 3 A £ ~ i T?T 2 2 5(T 3 )e- m ( T i +T2+T ^, 


where we have put D = 4 — 2e and A = T\T 2 + T 2 T 3 + T 3 T\. Also, the damping mass factor 
e -m 2 (T1+T2+T3) - g i nser teci in each Cj in the same way as in the one-loop calculation (Section |||). 
In the meantime, we shall introduce the notation e' = {e;e > 1} for C\ and Cq, in order not to 
mix it up with the (usual) infinitesimal parameter e > 0. This description is indispensable for the 
convergency of the C\ and Cq integrals, although we shall set e' = e after all, expecting the analytic 
continuation (see also Section 8-1-2 in |l6|]). The e' divergences are related to the divergences from 
the artificial mass term insertion, since C\ and Cq contain tadpole contributions, which vanish in 
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the m —* limit (in the sense of dimensional regularization) . 

In C5 and Cq, all the integrals are easy to perform, and hence we simply write down 



c 6 



(m 



2\-2e 



r(e + i)iy - 1) 



(5.4) 
(5.5) 



The rest of Cj are computed in detail in Appendix [C], and we only show the results as follows: 



Co 



c. 



+ 



+ 



m 



2\-2g 



B(e / -l,e + 2) + ( 



rn 



2\-2e 



1 

64 



r(e + 3)r(e-3) 



3-e 

4B(2 i) F(2,3 + e,^ i) - 3B(3 i) F(3,3 + e,^) 



2'4' 



2' 4' 



m 



T(2e + 1) 
(3 -e)(e- 2) 



B(e,i) 3 F2(l,2e + l,e;e-l,e + ^) 



^B(e + 1, ~) 3 F 2 (1, 2e + 1, e + 1; e - 1, e + |, i) 



m 2 )- 2e ^-r(3 + e)r(e - 3)B(3, -) F(3, 3 + e; -; ^) 



1, 



7 1, 



(34 



2'4' 
1 1. 



m2) - 2£4 - £ £M B (e, i) 3 F 2 (2e, 1, e; e - 2, e + |;^) , 

-C< + (^1^ + ^- 2) 1 } F (4 _ 2 + 3 _ 5 1 ; 
v ; 16 (e + 2)(e + 3) v 2 7 v ' 2 4' 

m2) - 2£4 - £r(2g) r(2-e) B( ^ 1 )3jF2(2) > 2£) £ . £ _ 1; e + 1 . i) 5 

^1^1^ 1 _ 3 _ 7 1 

; 32 (e + 2)(e + 3) v 2 J v ' ' 2 4' 

m2) - 2£ 1 4 - £r(2e) r|2^£| B(£ + ^ i )gF2(2j ^ e + 1; e _ ^ £ + 3 . 1 } 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



We can rewrite these expressions in terms of the hypergeometric function F = 2F1 only, and the 
results and their derivations are presented in Appendix |C] (see ( C.20 ), ( C.21| ), ( |C.38| ) and ( |C.39|) ). 
Since we could not find any convenient transformation formula from 3F2 to 2-P1 in the literature, 
we have established a transformation technique in Appendix [C]. (A more concise explanation can 
be found in Appendix |D].) 

Let us consider the expressions ( p.20|) , ( p.21| ), ( p. 3§| ) and ( p. 39 ). We now perform the Taylor 
expansions of the hypergeometric functions around e = 0, in order to see the divergence structures 
of the coefficients Cj. Here, we are only interested in 1/e terms in the sense of the MS scheme, and 
the hypergeometric functions which contribute to the desired pole terms are only generated from 
the following expansion: 

F(a, ae; 7 + be; -) = F(a, 0; 7; -) 



16 



+ aeF^ ( a> o; 7; I) + fcF^ ' 1 ) (a, 0; 7; i) + 0(e 2 ) , (5.10) 

where 

F^ l \a,f3; T ,z) d = d» df % F(a, # 7; z) . (5.11) 
These differential coefficients (for 7 7^ 0) are evaluated by 

F(a,0; 7 ;*) = l, F^fo 0; 7; z) = , (5.12) 

F^Xa&TiZ) = F^°\0,a;r,z) = -z 3 F 2 (l, 1, a + 1; 2, 7 + 1; z) . (5.13) 

7 

Here we again encounter the generalized hypergeometric function 3F2, however in the present 
case, it can be reduced to the ordinary hypergeometric function 2F1 through the following formula 
(derived in Appendix |D|): 

1 n 

3 F 2 (l,(3,n + l;2, r ,z) = -J2 F ( k >fcr,z), (n > 1, |z| < 1, 5R( 7 ) > 0). (5.14) 
Combining ( jTD| ) and ( p^) ; we have 

n 

i^V.OiTiz) = -£F(fc,l;7 + M , (n > 1, |*| < 1, »( 7 ) > 0), (5.15) 
T fc=i 

and thus 



F(n, ae; 7 + be; z) = 1 + e— J2 F ( k ^ h 7 + 1; z ) + °( £2 ) , (n>l, \z\ < 1, ^(7) > 0). (5.16) 

^ fc=i 

Owing to this formula, all coefficients in front of 1/e can be written in terms of 2-F1 and the gamma 
functions. After some algebra, we obtain 

5 + + , ^ = ^+(--^- + 0(1) , 

9 3 /e w 2 6e 2 V18 3 > e K ' 

-O(l), C ' e = ~ + 0(1), (5.17) 
where the overall factor (47r^ 2 ) 2e seen in (|5.l| ) and (|5^) is absorbed in Q; i.e., 

a = (4tt M 2 ) 2£ Q , (5.18) 

and we have defined 

?7i 2 

= 7e + hr - — ~ , 7 B = Euler const. (5.19) 
47T/ir 





1 






+ (- 


C 3 = 


1 




12e 2 


-f 1 

V72 




1 
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Therefore from Eqs. fl5,l| ), ( |5.2[ ) and ( 5.17| ), we obtain 



h[A] 
h{A] 



4.9, 



(4vr) 4 



+ 



11 

T 



45o 4 



(4vr) 



+ 



8/ ° m )l J d ° x o [-1^^) + ( 5 - 2 °) 
d D x f-^^ + O(e ) (5.21) 



and due to Eq. Q2.1| ) the renormalization part of the effective action (purely gluon parts) at the 
2nd order in T in our regularization is found to be 



r[.A] 



+ O(e ) 



(5.22) 






(a) 



(b) 



Figure 2: The "eight-figure" diagrams contained in the coefficients C5 and C6 



Finally, let us put a comment on what our results imply. Picking up C5 and Cq from Eqs. (|5.1j) 
and (5.2), and using ( |3.17| ), we extract the following quantities corresponding to self-energy parts: 

45 ab 



n 5 



ab 



9t 



(4tt) 4 
TT ab - ^ 



{k 8fu/ — kfj,k v )(— 6C5) , 
(A: 2 V - k^)(-WC' 6 ) . 



(5.23) 
(5.24) 



As briefly shown in Appendix [E], the H5 and 116 exactly coincide with the Feynman diagram results, 
if the coefficients of gluon kinetic terms are evaluated in the region close to the light cone k 2 — ► 0; 
in other words, if k 2 is much smaller than the mass parameter m 2 . (Note that Ca = 2 in the 
su(2) case.) Thus, it is very natural to expect that the other coefficients C±, C%, C3 and C4 should 
possess the same meaning. 



6 Conclusions and Discussions 



In this paper, we have explicitly calculated the gluon parts of the two-loop Euler-Heisenberg 
actions, which are organized at the level of an implicit formulation in the previous paper [11]. The 
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present results are still preliminary to reach a clear physical quantity such as /3-functions, however 
this paper is an important step toward a full-fledged extension of the world-line formalism to two- 
loop Yang-Mills theories. One of the main obstacles for this aim is the problem of how to integrate 
the proper time variables at higher loop calculations. Also, in the sense of field theory limit of 
string theory, this is an important problem: at the level of string theory, it is recognized as the 
problem how to perform the moduli integrals on a multiloop world-sheet. 

We have performed the path integrals in Section |], applying the world-line Green functions 
to all combinations of s«(2)-like charges (the k signatures; q.v. ( 4.11| ), ( (4.12| ), ( [4.14 ) and ( |4.15| )). 
Then extracting the parts corresponding to the wave function renormalization, we have been able 
to perform the world-line "moduli" integrals to reveal the divergence structure in Section |5[ In the 
standard method, it is difficult to perform the loop integrals containing a mass parameter, while 
in our case, all the integrations are carried out in Appendix |C], giving rise to (generalized) hyper- 
geometric functions as a result. This is certainly a significant point from a theoretical viewpoint, 
and hence we have mainly focused on the technical issue concerning the integrations. We should 
also note that the pseudo-abelian technique has worked out both at one- and two-loop levels, with 
reproducing the Feynman diagram results of Appendix [E]. We expect that our integration method 
will straightforwardly apply to higher order terms in T as well. 

Although we have verified that our results contain correct Feynman diagram contributions, the 
followings should further be investigated as a next step toward our goal: the coincidence between 
the present results and those by the Feynman diagram method can only be understood in the 
region close to the light cone k 2 — » (k 2 << m 2 ), where k^ are the external gluon momenta. 
On the other hand, as seen in Section ||, we do not have the restriction on k 2 at the one-loop 
level, in order to extract the /3-function. Similarly we shall encounter the same difference in the 
ghost loop calculations, and should clarify the reason for this kind of discrepancy. Related to this 
issue, another question is whether or not we can evaluate the pole structures of Cj for m = (or 
m 2 « k 2 ). As inferred from Eqs. ( |E.14j ) and ( |E.16 ), the integrals C{ might depend on the region 
of either k 2 « m 2 or not. However the present calculations do not indicate such a dependence, 
simply because Cj are not the Fourier modes of gluon two-point function. To clarify this point, 
one should compute the correlator of two gluon vertex operators (of bosonic field representation) 
along the outline of Appendix B in Ref. [O]; in this case, the formulation should be extended to the 
super world-line formalism in order to optimize the inclusion of the four point interactions involving 
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external legs. Anyway, in order to find the correct /3-function coefficient at two loops, we also have 
to add the contributions including the counter terms generated from one-loop divergences (in the 
massive formulation). 

Following ref. fl20|| a gauge symmetry breaking IR gluon mass m 2 (introduced there only as a 
device to seperate IR and UV divergences) requires a further counter term 

-Z x m 2 AlAl with z x = -^4~ (6.1) 

2 MA* * 167r 2 £ V ) 

and Ca = 2 in our case, cancelling m 2 dependent singularities in the MS dimensional regularization 
scheme. Insertion of this counter term 'gluon mass' into the one-loop contribution of order T 2 as 
found in Eq. ( 3.15| ), i.e. calculating with m 2 — > m 2 + Z x m 2 and expanding to first order in thus 



using 

e-W^(l-^r), (6.2, 
yields, e.g. a singular two-loop ^-contribution: 

°° 2 2 

-at = ^^(g-2) JdTT^e-^ C^^ T) Jd D x Tr L T 2 (6.3) 



2(4tt) d / 2 

„4^2 



|(i D xo f-ij>J>J + O(e ) • (6.4) 



(4tt) 4 

However in the background formalism there should be further counter terms including quantum 
as well as background fields. This has to be further analyzed in order to reproduce the usual 
/3-function coefficient. 

This paper concerns a theoretical interest, and is an important part of the ongoing long-term 
effort to find a way to getting over difficulties in the current calculation methods for higher loop 
amplitudes in field theory. Since the whole computation process in the world-line formalism looks 
completely different from the standard field theory calculations, the above questions are the mile- 
stones in the future study, and should be solved in order to make future practical applications 
successful. 
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A List of path reversal formulae 

This appendix is a brief note on the reversals of path ordering and proper time directions. The 
standard definition of the path ordering (the normal type) is 

J dTM T [x] I = E / dT l / • • • / dT « M n [a]'" M r n [*] . (A.l) 

and the anti-path ordering (used in [|ll|] for a certain reason) is 

J drM T N =£ jdnjdr 2 ---j dr n M Tn [x] ■ ■ ■ M T1 [x] . (A.2) 



Ta 



We here assume that the M r in the above two definitions are the same objects. 

The relations between the path and the anti-path ordering formulae are given by 

{rp v ae f T \ ea 

J drM T [x] \ = P* exp J J drM?[x] \ , (A.3) 
J ,w 1 ) vu, 



W^ pu [x;S,T a ,0]=Tv c 



X a P* exp J J drMj [x] 1 A e P* exp j ^ cfrMj [x] 1 , ( A.4) 

^ ) vp Ktc ) alx 

where the M T represents the transposition on both the color and Lorentz spaces. 

Note that the transposition makes an additional minus sign in front of the Ax term in M T . If 
we change the sign of gauge coupling g, this additional sign drops out, and the starting formulae 
presented in Sect. [2] follows from the previous paper directly. 

The other useful observation is related to path inversion: As can be seen from Eq. (^?j) by 
distinguishing r and r' carefully (see Eq. the following relation holds: 

(M T [x- l ]) a ^ = (M Ta+T ^ T [x]) b ^ , (A.5) 

where x~ l denotes the inverted path x, each of them defined by 

x ■ [r a ,Tp] -> R A , t\->x(t), (A.6) 

x' 1 : [r Q ,Tg] -> i? 4 , t i-> x' 1 (r) := x{r a + 773 - t) (A.7) 



respectively. As a consequence of Eq. ( A.5 ) we have 







T 




P exp < 


JdTM T [x] 


> = P exp < 


JdrM T [x- 1 ^ 











(A. 
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where again T means transposition in both color as well as Lorentz space. 

Describing path inversion by the x~ l symbol is a useful tool in performing world-line calculations 
(rather than changing the direction of r or other possibilities), because it nicely fits together with 
the following identities for world-line path integrals^]: 

x(S)=y 2 x(S)=yi 

J Vx F[x] = j Vx (A. 9) 

x(0)=yi ie(0)=j/2 

and consequently 

x(S)=y x(S)=y 

j>Vx F[x] = Jd D y J Vx F[x] = Jd D y J Vx F[x^} = j> Vx F^' 1 } , (A.10) 

x(0)=y x(0)=y 

where F is an arbitrary functional and the integrands -F[x _1 ] are to be understood as follows: for 
any path x integrated over within the path integral the corresponding path x _1 is to be constructed 
(in thoughts) and the functional is to be evaluated for this path Thus, in the path integral, 
depends on x. Note that we also have 

/ [Vx] s F[x] = / Vx e o F[x] = Vx e 

x(0)=yi x(0)=yi x(0)=y 2 

x(S)=yi f, 12 x(S)=yi 
f ~ J dT 4 X f 

/ Vx e o F[x~ l ] = / [Vx] s Fix" 1 ] , (A.ll) 

x(0)=y2 x(0)=y2 

i.e. our usual bracket notation for the free path integral part is not affected. 

Using these relations and identities one can easily verify that the gluon propagator in a back- 
gound, given by the world-line representation pi 

oo x(T)=X! t T \ ab 

A^( Xl ,x 2 ) = J dT J [Vx] T Pexp J JdrM T [x] 1 , (A.12) 

x(0)=x 2 ^0 J ^ 

satisfies the following property: 

A%(x 1 ,x 2 ) = A%(x 2 ,x 1 ) . (A.13) 



2 These identities can easily be verified from the path integral discretization. 



22 



B The derivation of I\[A] 



In this Appendix, we show some details for the computation in Section 4.1. First, we need 
various values of the quantity ( 4.17 ) at (r, r') = (0, T3) for a,b = 1,3, and those are given by the 
following. The necessary derivatives of the Green functions are 

/ 2T 2 T 2 e^-Ti) \ 



d T d T ,gll(T,T';+,0,-) 



drd T ^(T,r';0,-,+; 



d T d T ,gli(T,T';+,-,o) 



^^3(r,r';+,0,-) 



t=o,t'=t 3 y sin TT X sin TT Z + TT 2 sin .F(Ti + T 3 ) 

_ / 2T -sm^e-^ 

T= o,r'=T 3 ~ ^ sin TT 2 sin TT Z + J^Ti sin T{T 2 + T 3 ) 

/ 2J"sinJ"T 2 e-^ Tl 
r =o,r'=T 3 ~ ^ sin J=T X sin ^T 2 + ^ sin + T 2 ) 

/ 2T 2 T 2 e i2TT * sin J"Ti cosec ^T 3 
t=o,t'=t 3 ~ I sin J^Ti sin^"r 3 + TT 2 sin^"(r 1 + T 3 ) 



fJ,V 



[IV 



(B.l) 
(B.2) 
(B.3) 



+ 2 1 L 6(T< 



£U> t ^(t,t';0,-,+) 



sinJT 3 , ^ 

__ / 2riJ" 2 e~ i2 - FT3 sinJT 2 cosecJ c "T3 
T=0 ,r'=T 3 ~~ I sin ^T 2 sin J^T 3 + ~TT 1 sin .F(T 2 + T 3 ) 



(B.4) 



+ 2 1 L «S(T 3 



sin J-T 3 



(B.5) 



cU^(t,t';+,-,0) 



2To 1 sin TT\ sin J^"T 2 



T =o,r'=r 3 I sin TT X sin J^T 2 + TT 3 sin .F(Ti + T 2 ) 



The normalization constants are given by 

^(+,-,0) = (47r )-D det l/2 7 



+ 2(5(r 3 )--)(i L V . 



^ 2 



(B.6) 



|d D x , (B.7) 



sin J^Ti sin TT 2 + .FT3 sin F(T\ + T 2 ) 

and the other values can be obtained by exchanging the labels a on K a and T a simultaneously; for 
example, 



^(+,0,-) = ^(+,-,0) 



Note also 



AA(+'+'°) =J\f { ~~' 0) =M {± ' T ' 0) . 



(B. 



(B.9) 



Now, plugging in ( |B.1| - |B.6| ) and the corresponding normalizations, for each line of Eqs. (4.14) and 
( 4.15| ) we get an expression of the following structure: 



series in J 7 ) + Tr^ (power series in J 7 ) ■ Tr^ (power series in J 7 )] . (B.10) 
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Because of the antisymmetry of J 7 , only even powers of T contribute to the traces. Furthermore 
we have the property 

and thus ( |B.10| ) is invariant with respect to T — ► —J 7 . Therefore the .F — ► — F terms in Eqs. ( 4.14| ) 
and ( [4.15| ) just give a factor of two. 

Using this fact, the expressions (B.1- [B.6| ) and the corresponding normalizations we obtain from 
Eqs. ( 4.14 ) and ( 4.15| ) straightforwardly: 



and 



with 



(4tt)- d jdT l dT 2 dT z det 1 / 2 



J 7 ' 2 



A 



(2) 
T ■ 



F 2 T 2 r 



A 



(2) 
T 



-2 cos T(Ti - T 3 ) cos 2T{T X + T 3 ) + cos JF(T 3 - T{) 



+ cosF(Ti - T 3 )Tr L cos 2T(T X + T 3 ) 



+ 2Tr L 



A 



(2) 
T 



2 sin FT 3 sin 2T(T 1 + T 3 ) - cos F(2Ti - T 3 



+ cos F(2Ti + T 3 )Tr L cos 2FT 3 



ri[A] = 



X - (4vr)- D [dTtdndTs det 1 / 2 



T 2 



A 



(2) 



sinFT 3 

+ 2 cos FT 3 Tri cos 2.F(7i + T 3 ) 
sin TT\ sin FT 3 1 



o v ~ ^ ' * 

-2 cos F(2Ti + 3T 3 ) + cos T{2T X + T 3 )Tr L cos 2JFT 3 



Tr L ( 



+ Tr L 



a<?>t 2 

F 2 r 2 sinFri 
A^ } sinFT 3 



T 2 ' 



cos 2F(T 3 - Ti) - cos 2FT 3 Tr L cos 2FTi 



2 cos 2F(Ti + 2T 3 ) - cos 2T{T X + T 3 )Tr L cos 2FT 3 



cos 2FT 3 Tr L cos 2^"(Ti + T 3 ) + 5(T 3 ) 2(1 - D) Tr L cos 2FIi 



+ <J(T 2 ) Tr L cos 2F(Ti - T 3 ) - Tr L (cos 2TT X ) ■ Tr L (cos 2FT 3 ) 



(B.12) 



(B.13) 



A {2) = sin FTi sin JFT 3 + FT 2 sin JF(Ti + T 3 ) . (B. 14) 

The sum of (B.12) and (B.13) gives Ii[A] by definition, and is certainly equivalent to Eq. (4.21). 
In order to obtain Eq. ( [4.21 ) itself, we should further take the following modification into account. 
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Exchanging T 2 and T3, and using the relations followed from ( |B . 14] ) 



T 



sin FT 3 & 



sinJTTi + 



^rasin^Ti + T 3 
AP smFT 3 



(B.15) 



1 sin TT\ sin TT 3 
T 2 aJ¥ 2 



A 



sin.F(Ti + T 3 



(B.16) 



we finally arrive at the full expression shown in Eq. ( [4 . 2 1| ) . 



C Computational details of C{ 

In this Appendix, we show the details of how to perform all the integrals in C\, C 2 , C3 and C4. 
Let us first perform the T3 integrals in Cj. Applying the following formula to the T3 parts in (15. 



a - pt {l + at)- v dt = tf- x oT v e p l a T(l - u; p/a) , (a > 0) 



and then transforming the T\ and T 2 integrals with 

/•oo /"OO rl 

/ dT 1 dT 2 f(T 1 ,T 2 ) = dTT / duf(T(l — u),Tu) , 

JO JO j0 

we obtain the followings for z = 1 and 2 (further using m 2 T = t) 



00 /■! 



Q = (m z )-' Z£ I dt / dufi(u)t 2+£ e- t+tu(1 - u) T{e - 3; tu{\ - u)), (i = 1, 2) 
j j 



and for i = 3 and 4 

Ci = - 



00 /•! 



■(m ) / dt / dufi(u)t e e 
Jo Jo 



£ -t+tu(l — U) 



(3-e)tr(e-3;tu(l-«)) 



- «)r(e - 3; tu(l - u)) + t 2 u(l - u)—T(e - 3; £) 



where we have defined 



and 



f = tu(l - u) 



i-1 



G 



(C.l) 



(C.2) 



(C.3) 



(C.4) 



(C.5) 



(C.6) 



f l (u)=u*- 2a (l-u) a+5 - i , a 

with Gauss' integer symbol [ ]q. In the following, we evaluate ( |C.3j ) and ( C.4 ) separately because 
we shall proceed on different technique and formulae. 



25 



C.l d and C 2 



We now consider the t integration in (C.3). In the first place, it can be integrated in terms of 
the formula (6.455.1 in pi): 



(C.7) 



t ^e-^Y{v, at) dt = "^t+l F 0-> H + v, fj, + 1; -^-), 
o ji[a+p)^ +u a + p 

Re (a + p), Re fi, Re (fj, + v) > . 



Then applying the formula (9.131.2 in |i9| ), 

F(a,(3; T ,z) = ^^J^(l - z)^^ F(j - a, 7 - /3; 7 - a - + 1; 1 - z) 

+ r(7)r( 7 - « - P) F ^ (3 . a + f3 _ 7 + 1 . 1 _ z)i (a8) 



r( 7 -a)r( 7 -/3) 



we have 



d = (m 2 y 2£ f 1 dufi{u)Y[e + 3)r(e - 3)F(4 - e, 3 + e; 4 - e; «(1 - u)) 
Jo 

-ar(2e) 



+ (m 



2\-2e 



dufi(u) u(l — it) 



3-e 



F(2e,l;e- 2; «(!-«)) . 



(C.9) 



In order to perform the u integrations, we consider the following change of variable: cutting the 
integration region (0, 1) in half, define u + (it-) for larger (smaller) values of u, and apply 



u± = -(l±VT=i) • 

With this change of variable, the following formula holds for an arbitrary function H{u) 

1 If 1 dv 

dufi(u)H(u(l -u)) = - —^= gi (y)H(y/A) , 

o 4 Jo v 1 - y 

where we have defined gi(y) for all i: 

9i(y) = /i( 

and we hence have 



l + ^r^ ,1-yr^ 



(CIO) 



(C.ll) 



(C.12) 



9i(y) = T^y( 4 - 3y), ff 2 (y) = T^? 
lb lb 

53(y) = o2/(2 - y), g±{y) = ^y 2 . 



(C.13) 
(C.14) 



At a glance, one may realize that the y integration can be performed by the formula (7.512.12 [|l9| ]) 



(1 - t^'H^Fia, (3; 7; zt) dt = B(/i, v) 3 F 2 (a, (3, + v\ z), 



(Re [i, Rev > 0, \z\ < 1) 



(C.15) 



2b 



and this is exactly the way we obtain the generalized hypergeometric 3P2 expressions ( |5.6| ) and 
& 

In order to further derive hypergeometric 2-P1 expressions, we rather notice the special case 
7 = v in ( C.15|) . Before applying (|C.15|) , we adjust auxiliary variables of 2 Pi until the special case 
is applicable, with using the formula (9.137.18 (|Tg[|) 



F(a, 0; 7; 2) = l-^F{a, (3; 7 + l;z) + -F(a + 1, ft 7 + 1; z) 

7 7 



(C.16) 



Before writing down the final results for C\ and C2, we here remark on the way how C\ contains 
the e' dependence. Using the following formula (9.137.11 ]l^]) with (3 = 



7-F(a, (3; 7; z) - -yF(a, (3 + 1; 7; a) + azF(a + 1, /? + 1; 7 + 1; 2) = 



(C.17) 



and then decomposing the second term in (|C.9|) 



F(2s, 1; e - 2; |) = 1 + | F(2e + 1, 1; e - 1; |) 

we extract the following integral from C\: 

d ^ 4" £ r(2e) /' 1 ^- 2 (4-3y)(l-y)- 1 / 2 ( iy 
Jo 



(C.18) 



/1 



r(2e) / X u £ ~ 2 (l - u) £+1 du = r(2e)B(e' - 1, e + 2) 
Jo 



(C.19) 



Evaluating the l\ contribution separately, we obtain the results written in terms of 2 -Pi only: 
Ci 



(™ 2 )- 2£ ?^B(e' - 1, e + 2) + (m 2 )- 2£ -Lr(e + 3)r(e - 3) 
3 — e o4 



4B(2, 1) F(2, 3 + e; 1 1) - 3B(3, ±) F(3, 3 + e; I; 1) 



+ (m^ £ 4 



2,-2e A -e r(& + 1)B( £ , ±) 



(3- e )(e-2)(e-l) 



(e - 2)F(1, 2e + 1; e + i; i) + F(2, 2 £ + 1; e + 1; ~) 



[m 



2\-2e A -e 



4" e - 



3r(2e + l)B(e + l,i; 



4 (3-e)(e-2) 



e-2 



F(l,2e + l;e + ^;i) 



(C.20) 



Co 



1 



1 



7 1, 



(m } 64 r(3 + £)r(£ " 3)B(3 ' 2 )F(3 ' 3 + £; 2 ; 4 } 



1-e 

2 11, 
F(2,2e; £ + -;-) + 



F(l,2e;e + 



2' 4' 



2-e 



2 11, 

2 '4 / ' (2-e)(3-e) F ^ 2e;£+ 2 ! 4' 



(C.21) 
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C.2 C 3 and C 4 

In this case, we split Q; i = 3,4 as 



where 



Ri 
Rs 



Ci = -(m 2 )- 2£ (3 - e)Ri + i? 2 + R3 , 



1 /-oo 

fi(u) / t 1+e e-' + ' u(1 - u) r(e - 3; tu(l - u)) dt , 
Jo 

1 /-co 

du fi(u)u(l - u) / t 2+£ e- t+tu(1 - u) r(e - 3; tu(l - u)) dt 
Jo 

1 roo f) 

dufi(u)u{l-u) / t 2+E e- t+tu ^- u ^T(e-3;^)dt . 
Jo o<i 



(C.22) 



(C.23) 
(C.24) 
(C.25) 



First, we integrate R2 using ( |C.7|) in the same way as done in (|C.3[) : 

R 2 = / fi(u) \ «(1 - u) 1 6 ~ 2e; e + 4; 1 - u(l - «)) du. 



e + 3 Jo 

Second, we notice that the incomplete gamma function is related to the Whittaker function: 



(C.26) 



r(i/,z) = z^e-'^W^i u{z), 



(C.27) 



and W K;U (z) satisfy 



W K Jz) = W.-Jz) 



Z 



zd z W^(z) = (| - ^W^z) - W K+1 ^(z), 
W K ^(z) = z 1/2 W K _ 1/2jfM+1/2 (z) + (_-«- //)W K _i )M (z) 



(C.28) 
(C.29) 

(C.30) 



Applying the formula ( C.27|) and its derivative to R\ and R3, one can prove the relation 



R3 — — g — ^ ~ ^ 4 



(C.31) 



where 

i? 4 = / d«/<(«) 2 / t 3£/ V* + * u(1 -" )/2 d £ WV4 ^s(e)dt . (C.32) 

Jo L J Jo 2 ' 2 

We here remove the derivative d^W from the r.h.s. of Eq.(C.32), making use of the relation ( C,29| ): 

1 £ _4 fl r ,£^4 /-OO 



-TLA — XL9 — 

2 2 



dufi(u) u(l — u) 



t _ e _ le _ t+tu( i_ u)/ 2 w ^ ^3(£)cft . (C.33) 
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Adjusting the indices on W on the r.h.s. of flC,33j ) in terms of the recursion relation ( C.30| ) (with 
k = (e — 2)/2, fj, = (e — 3)/2), we then apply the formula 



t w W K> ^(at)e- pt dt 



r(u + fi + |)r(i/ |)a^ +1 /2 

r(i/ - k + 2){p + a/2) v+ ' 1+3 / 2 

3 1 

xF^ + /i+-,/j-K+-;i/-K + 2; 



(Rez/ + - - |Re/i| > 0) 



2p — a 
2p + aV' 



(C.34) 



and the r.h.s. of ( C.33| ) becomes 



2 2 v ' Tfe + 3 













u(l—u) 


/o 





e-2 



Substituting this expression for in Eq.( |C,3i] ), which should then be inserted into Eq.(C.22), and 
thereby eliminating R±, we calculate Cj as follows: 

r h 6 _2_,. , r -1 



C; 



= {m 2 )- 2£ du f l {u) u{l - u) r(2e) 

JO L J 

+ ^— F(l,2e;e + 3;l-«(l-«)) 
(m 2 )- 2e r(2e) /•! 



3 + e 



F(l,2e;e + 4;l-u(l-u)) 



, -y)"^(-) e " 2 F(2,2e;£ + 4;1- hdy 

4(2 +e)(3 + e) Jo 4 ^ 



(C.36) 



where we have used the transformations ( C.ll| ) and ( C.16| ) at the second equality. We now transform 
the argument of the hypergeometric function from 1 — y/4toy/4 through the formula ( |C.8| ), and 



thus giving rise to F(4 — e, 2 + e; 3 — e; y/4) and F(2, 2e; e — 1; y/4). (If we apply ( C.15 ) at this 
stage, we obtain the generalized hypergeometric 3F2 expressions ( |5.8| ) and (p\9|).) 

For the purpose to reduce 3F2 to 2F1, we apply ( C.16Q twice to the latter F (two of them 
explained right above), and apply the following formula to the former F (setting a = 7 = 3 — e, (3 = 
2 + e): 

R n, - R 

(C.37) 



F(a + 1, 0- 7; z) = -F(a, R + l- T)Z ) + ?—?-F(a, 0; 7; z) . 
a a 



After that, we are able to integrate them by using ( |C.15| ) (with 7 = 1/), and the results are therefore 



C 3 



(m 2 ) 



2\-2e 



1 r(4 + e)r(g-2) 
16 (e + 2)(e + 3)(3-e)~ v "'2 



B(2,-l 



(2 + E )F(2,3 + e;|i) 



+(l-2e)F(2,2 + e;|i; 

+ (m 2 )-^r(2 £ ) r(2 - £) B( ^ e) 

+2F(3,2e;e + i;i)] - C 4 , 



(e-3)F(2,2e;e + 



1.1, 

2' 4' 



(C.38) 
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Ca 



2 2£ 1 r(4 + e)r( g -2) 1 r 7.1, 

(m } 32( £ + 2)( £ + 3)(3- e ) B(3 ' 2 H (2 + £)F(3 ' 3 + £ ' 2 ' 4' 

+(l-2e)F(3,2 + e ;I;i) 



+ (mT r(2 £ ) 



r(2-e)B(i,e + l) 



r(4-e) 2e(e-l)4 e 



( £ _ 3 )( £ _2)F(2,2 £ ; £ + ^;i) 



+4(e - 3)F(3, 2e; e + | i) + 6F(4, 2e; e + |; 1) 



(C.39) 



D The proof of Eq.(5.14) 



We show the outline of how to prove the formula ( |5.14D . Basically we follow the same technique 
as we used in Appendix [C], with using ( C.15 ) and ( |C,16|) for the aim of rearranging 3F2 into 2 -Pi- 
We evaluate the l.h.s. of Eq. ( p-15| ) in two ways: one is the direct result (the r.h.s. of the formula), 
and the other is a combination with ( |C.16 ). 

On the one hand, setting 7 = v — n in ( C.15| ), we have 
-1 



(D.l) 



(1 - ty L t u ~ L F(a, 0;u-n; zt) dt = B(/i, v) 3 F 2 (a, (3,v;v - n, \x + v\ z), 
(Reu,Reis > 0, \z\ < 1) . 



On the other hand, using ( C.16| ) n > 1 times, the integrand in Eq. ( D.l| ) can be expanded as 



v — a — k ^-^ -r-r a + p — I n 



En 



v — k v — a — p\r 

k=l r=0jj=l 1 \ / 



F(a + r,f3;v;zt) , 



(D.2) 



where we define ]T2=i = 1- Then applying ( C.15|) (with the case 7 = v) to this expression, the 



l.h.s. of (D.l) becomes 

n 

B^n 

k= 

We thus have the equality 



1 n r , i / \ 

v — a — k ^-^ rr ol + p — 1 I n 



En 



, , v — k ^. J - J \v — a — p\r 

k=l r=0p=l 1 \ / 



F(a + r, f3; (J, + v; z) 



v — k '—l ■*•■*; v — a— p\r t 

k=l r=0p=l r \ J 



Putting a = 1, v = n + 2, and using the identity 



p In 



n — - 

i n + l — p\r 

p=i ^ \ / 



1 - 



we therefore have proven the formula ( ^.14[) : 



(D.3) 



„ , . -A- v — a — k n «+» - 1 fn, 
3 F 2 (a,/3,i/;^-n,/x + i/;z) = jj r ^l]: — r — 1 1 _ ) ^(a + r, /?; /i + z^;z) . (D.4) 



(D.5) 



? F 2 (l,/3,n + 2;2, / u + ^;z) = — — V F(r + 1, /3;» + 1/ 

n + 1 e -*L 



r=0 



(D.6) 
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which holds for KQu) , 3ft(V) > 0, \z\ < 1 and n > (we have proven it for n > 1, but the n = case 
is trivial). 



E Feynman diagram results 

In this Appendix, we present some Feynman diagram calculations in reference to the results 
obtained by our method. We follow the same notations as Refs. [17, 15] in the Minkowski space, 
however use the massive propagator in the (background) Feynman gauge: 

a b 5 ab g^ 
A A A A A A A = (-i) . 

V V V V \J V ,2 2 

We only deal with the parts which contain the "eight-figure" vacuum diagram. When the 
massive propagator is introduced, the tadpole contributions remain (see the diagrams (a) and (b) 
in Figure 2). After some calculations, the tadpole diagram (a) reads 



n (a) - a 



/~<2 rafe 

,4 

( 47r )2-e ' 



(m 2 ) 1 - £ (3-2e)r(e-l) 8(k 2 g^ - k^k u ) J x + D{k^k v J x + 4k u J 2 + 4J 3 ) , (E.l) 



where 



■h 

J2 



d D p 




1 


(2ir) D i 


{(p + k) 2 


— m 2 ){p 2 — m 2 ) 2 








(27T) D i 


((p + k) 2 


— m 2 ){p 2 — m 2 ) 2 


d D p 




VlxPv 


(2ir) D i 


((p + k) 2 


— m 2 ){p 2 — to 2 ) 2 



(E.2) 
(E.3) 
(E.4) 



and these are equal to 



T r ( g+1 ) • 



(4vr 



/ r ( £+1 ) . • 



(4vr 



■h 



r(e + i) , i . 1 7 7 • \ 

(47r)2-e ^2^^ j3 2 » J 2 ' ' 



(E.5) 



with 



Ji 

32 

h 



(2 2 / 2\\"^ 1 " 

m — k (x — x )) dx, 

l . 
x(l — x)(m 2 — k 2 (x — x 2 ) 

i , 

x(m 2 — k 2 (x — x 2 ) 



1 "■ ^ dx, 



£ = \(m 2 )- £ -\ek 2 (2j 2 -j l ) 



(E.6) 
(E.7) 
(E.8) 



31 



Thus Eq,( |E.l| ) becomes 



n 



(a) 



ab 



ri2 xab 

g 4 A .. (m 2 ^- e 



(3 - 2:)1'(: - + 1 ) j(m 



( 47r )4-2e 

+{k 2 9^ ~ k^K)((D - 8)ji - 2Dj 2 , 



2\-£ 



and the other tadpole diagram (b) is calculated as 



/~<2 sab 
4 U A° 



(m 2 ) L - 2£ gilu (3 - 2e)DT{e - l)r(e) . 



As a result, the sum of these tadpole contributions takes the transversal form: 



ab 



rrl a)C +nv 

/~i2 xab 
A ^ A° i 2\1— f 

g , . ,a n_ (m u £ 



r(6) ( 



(3 - 2e)r(e - l)r(e + l) 



(4tt) 4 - 2£ 

x - fc^) f (£> - 8)ji - 2Dj 2 



When k 2 — > (|/c 2 | << m 2 ), the quantities j$ behave as 



2\-l-e 



J2 = ^(- 2 )- 1 - £ 



and when m 2 — ► 0, they reduce to 



h = ~ 2 (-k 



2\-l-e 



B(-e, -e) , j 2 = (-A; 2 )- i - £ B(l - e, 1 - e) , 



j3 = -(-kT £ m-e,l-e) . 



Therefore in these limits Eq. flE.li] ) behaves 



{ C 2 A 5 ab 



HtJ = glj^r(k 2 g^- k^){ e 



— (k 2 ^0) 



(m = 0) 



+ O(l) 



Next, the diagram (c) amounts to 



C 2 A 5 ab 



n (e) ^ = Sgo^yr(k 2 g^ - ^^)exp[2 e ln(47r M 2 )]r 2 (e)(2j 3 ) 2 



where j$ is given by flE.8p . Substituting Eqs. ( E.12 ) and ( |E.13| ) for this js in each limit, 
the following limits of ( |E.15| ): 



6 12 

o H Pm 



(fc 2 0) 



(47T) 



6 24 12 
--^ 1 p (m = 0) 



+ 0(1) 



e e 
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where p m is given by ( |5.19 ), and 

P = 7E + ln^. (E.17) 

We therefore conclude that our results coincide with the Feynman diagram calculations in the 
situation k 2 -> (m / 0) (q.v. flpg) and (U|)): 



xab 

= 9 4 07 ^(k 2 g , u - fc M fc„)(-10C£) = n 6 J , (E.18) 



(4vr) 4 



n (e) ;, = 5 4 t£u ( fc2 ^ - = DbJ • (E.19) 
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